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THEORY  OF  VHF  SCATTERING  BY  FIELD-ALIGNED 
IRREGULARITIES  IN  THE  IONOSPHERE 


1.  INTRODUCTION 

Long  range  propagation  of  Very  High  Frequency  (VHF)  radio¬ 
waves  via  scattering  from  ionospheric  irregularities  is  of  impor¬ 
tance  to  the  design  of  Meteor  Burst  Communications  Systems  for 
the  following  reasons:  i)  it  can  be  used  as  a  means  to  communi¬ 
cate  at  distances  of  a  few  hundred  kilometers  in  the  absence  of 
meteor  trails  provided  the  communications  waveform  can  cope  with 
the  'distortions'  introduced  by  the  propagation  medium**  and  ii) 
when  present  along  with  meteor  trails,  it  can  interfere  with  the 
meteor  scatter  signal.  In  either  case,  it  is  important  to  know 
the  conditions  under  which  scattering  from  ionospheric  irregular¬ 
ities  will  occur,  the  strength  of  the  scattered  signal  and  other 
parameters  such  as  the  coherence  bandwidth  or  delay  spread  and 
the  coherence  time  or  Doppler  spread. 

Scattering  from  ionospheric  irregularities  can  in  principle 
occur  almost  anywhere  on  earth.  However,  it  has  been  most  fre¬ 
quently  observed  at  high  latitudes  (Collins  and  Forsyth,  1959] 
and  near  the  equator  (Ferguson,  1981]  while  very  weak  signals 
have  been  reported  at  mid-latitudes  [Bailey  et  al.,  1955].  The 
reason  for  this  is  that  the  scattering  of  VHF  radiowaves  at  high 
latitudes  is  associated  with  the  radio-aurora  which  produces  ir¬ 
regular  ionization  at  E-layer  heights  in  the  ionosphere.  The 
scattering  phenomenon  is  known  as  auroral  scatter.  The  equato¬ 
rial  observations  are  associated  with  nightime  field-aligned  ir¬ 
regularities  in  the  ionization  of  the  F-region  of  the  iono¬ 
sphere.  Because  the  phenomenology  of  the  high-latitude  radio¬ 
aurora  emulates  that  of  a  nuclear  burst  disturbed  ionosphere, 
strong  scattering  from  field-aligned  irregularities  at  VHF  can 
also  be  expected  to  occur  at  lower  latitudes  after  a  high- 
altitude  (100-500  km)  nuclear  burst. 
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In  Section  2  of  this  report  we  develop  a  theory  for  VHF 
scattering  by  field-aligned  ionospheric  irregularities  which  de¬ 
scribes  the  conditions  under  which  such  phenomena  occurs. 
Results  which  illustrate  the  geometrical  (aspect)  sensitivity  of 
the  scattering  by  field-aligned  irregularities  at  E-layer  heights 
(90-140  km)  in  the  auroral  region  are  presented  in  Section  3. 
The  theory  developed  is  very  general,  though,  and  also  applies  to 
scattering  by  field-aligned  irregularities  induced  by  high- 
altitude  nuclear  bursts  at  arbitrary  latitudes.  The  main  differ¬ 
ence  between  auroral  scatter  and  scattering  at  lower  latitudes  is 
that  at  high  latitudes  only  energy  scattered  by  irregularities  in 
the  E-layer  (90-140  km  heights)  can  be  received  by  a  terminal  on 
the  surface  of  the  earth.  Energy  scattered  by  irregularities  in 
the  ionization  of  the  F-layer  (200-500  km)  does  not  intersect  the 
surface  of  the  earth  because  of  the  directionality  caused  by  the 
anisotropy  of  field-aligned  irregularities.  At  lower  latitudes 
field-aligned  irregularities  in  the  E  and  F-regions  contribute  to 
the  total  received  scattered  field.  However,  the  scattering 
cross-section  of  the  irregularities  in  the  F-layer  is  greater 
than  that  of  the  E-layer  irregularities  so  that,  at  lower  lati¬ 
tudes,  most  of  the  received  energy  is  due  to  F-layer  irregulari¬ 
ties  . 

In  addition  to  the  aspect  sensitivity  (geometry)  of  field- 
aligned  scatterer  we  also  present  in  Section  3  results  which 
illustrate  polarization  effects,  antenna  height,  pointing  and 
beamwidth  effects,  and  frequency  dependence.  The  theory  needed 
to  determine  polarization  effects  is  discussed  in  Appendix  A. 

2.  THEORY  OF  FIELD-ALIGNED  SCATTER 

The  scattering  of  VHF  radiowaves  by  field-aligned  irregu¬ 
larities  of  ionization  in  the  ionosphere  has  been  treated  in  the 
literature  as  specular  scatter  (Stathacopoulos  and  Barry,  1974] 
or  as  scatter  from  long  cylindrical  irregularities  whose  ratio  of 
length  to  diameter  is  large  (Booker,  1956]  . 
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The  specular  scatter  theory  predicts  the  reception  oE  a 
scattered  signal  only  in  those  regions  where  the  specular  condi¬ 
tion  (to  be  defined  later)  is  satisfied.  This  restricts  the 
areas  of  the  world  in  which  field-aligned  scatter  would  be  ob¬ 
served  for  a  given  scattering  volume.  For  example  according  to 
the  specular  scatter  theory,  scattering  from  irregularities  at  an 
altitude  of  250  km  or  higher  (F-layer  scattering)  should  only  be 
observed  when  the  ground  based  transmitters  and  receivers  are  lo¬ 
cated  at  geomagnetic  latitudes  below  40®.  Similarly  scattering 
from  irregularities  at  altitudes  less  than  120  km  (E-layer  ir¬ 
regularities)  should  occur  only  when  the  transmitter  and  receiver 
are  at  geomagnetic  latitudes  below  62°.  However  there  is  experi¬ 
mental  evidence  that  shows  that  scattering  from  the  radio-aurora 
can  be  observed  in  the  neighborhood  of  the  auroral  zone  (geo¬ 
magnetic  latitudes  between  67®  and  73®  under  undisturbed  condi¬ 
tions),  and  sometimes  even  north  of  it  where  the  specular  condi¬ 
tion  cannot  be  satisfied  (Collins  and  Forsyth,  1959).  On  the 
other  hand,  there  is  also  evidence  that  the  scattered  signal  is 
stronger  and  can  be  observed  more  often  in  areas  where  the  specu¬ 
lar  condition  is  approached.  This  type  of  behavior  is  predicted 
when  the  irregularities  are  modeled  as  cylinders  whose  longitu¬ 
dinal  dimension  is  at  least  four  times  greater  than  their  dia¬ 
meter  and  such  that  their  axis  are  aligned  with  the  earth's  mag¬ 
netic  field.  This  is  essentially  the  theory  of  anisotropic  scat¬ 
tering  proposed  by  Booker  [1956]  to  explain  observations  of  scat¬ 
tering  from  the  radio-aurora. 

In  this  section  we  will  develop  a  model  for  calculating  and 
predicting  the  path  loss  of  anisotropic  scatter  given  the  geogra¬ 
phic  locations  of  the  transmitter  and  receiver,  and  the  height 
and  geographic  region  of  the  scattering  volume  where  the  field- 
aligned  irregularities  are  contained.  Since  the  scatter  from 
anisotropic  irregularities  is  aspect  sensitive,  it  is  always  use¬ 
ful  to  know  the  receiving  locations  where  the  specular  scatter 
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condition  is  satisfied.  Therefore  we  begin  with  a  derivation  of 
the  specular  scatter  condition. 


2.1  SPECULAR  SCATTER  CONDITIONS 

In  many  instances  it  may  be  necessary  to  determine  loca¬ 
tions  where  specular  scatter  from  field  aligned  irregularities 
may  be  expected  given  the  location  of  the  transmitter  and  the 
scatterer . 

The  condition  for  specular  scatter  from  elongated  irregu¬ 
larities  is  that  the  component  of  the  scattered  wave  vector, 

a 

Ug,  parallel  to  the  longitudinal  axis  of  the  irregularity, 
uf,  be  equal  to  the  component  of  incident  wave  vector, 
u^,  parallel  to  the  longitudinal  axis  of  the  irregularity,  i.e., 


u 


0 


u. 


(1) 


A  A  O  A  A 

where  a  •  b  is  the  dot  product'  between  vectors  a  and  b,  and 


u.  *  (P  -  T)/R  »  unit  vector  from  transmitter  at  T  to 
1  c  scatterer  at  P. 

A 

Uq  *  (R  -  P)/R  *  unit  vector  from  scatterer  at  p  to 

r  receiver  at  R. 

A 

uf  *  unit  vector  parallel  to  earth’s  magnetic  field  at  the 
scattering  point  P. 


The  incident  wave  vector  u.  can  be  decomposed  into  a  com- 

A  X  A 

ponent  parallel  to  u^  and  perpendicular  to  u^ 


as 


»  cos  a  Ug  +  sin  o  uN 


(2) 
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where 


uQ  «  coso  uf  +  sina  cos0  uN  +  sina  sin0  (u^x  uN )  (4) 


where  -*  <  0  <  *,  i.e.,  specular  scatter  due  to  an  elongated  ir¬ 
regularity  is  possible  within  a  cone  defined  by  the  angles  o  and 
3.  Figure  1  shows  two  typical  scattered  rays  which  intersect  the 
earth  at  possible  receiver  sites  Rj  and  R2«  The  locus  of  all 
such  intersection  points  defines  a  curve  on  the  surface  of  the 
earth  where  strong  scattered  signals  can  be  received.  Field- 
aligned  scatterers  are,  of  course,  not  concentrated  at  a  point, 
as  implied  in  Figure  1;  however  the  coverage  provided  by  many 
distributed  scatterers  may  be  calculated  by  summing  the  contribu¬ 
tions  from  each.  We  will  now  describe  how  to  find  the  locus  of 
scattered  rays  which  intersect  the  surface  of  the  earth  given  the 
geomagnetic  latitude,  0t,  and  geomagnetic  longitude,  of  the 

transmitter,  and  the  altitude,  h,  geomagnetic  latitude  0,  and 
geomagnetic  longitude,  $,  of  the  scatterer. 

Define  a  cartesian  coordinate  system  such  that  the  trans¬ 
mitter  location  is  given  by  T  =  a(cos0tcos$t,cos0tsin$t,sin0fc) 
and  the  scatterer  location  by  P  =  r(cos0  cos$,  cos0  sin$,  sin0) 
where  a  is  the  radius  of  the  earth  and  r  *  a  +  h.  We  want  to 
determine  the  geomagnetic  latitude,  0r,  and  longitude,  $r,  of  all 
locations  R  which  satisfy  the  specular  scatter  conditions  ( 1 )  — ( 4 ) 
such  that  | R j  *  a,  i.e.,  we  want  to  find  all  specular  scattered 
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rays  which  intersect  the  surface  of  the  earth.  We  will  assume 
that  the  incident  and  scattered  rays  propagate  in  straight  lines 
so  that  the  results  are  only  applicable  at  VHF  where  ray  bending 
in  the  E-layer  of  the  ionosphere  is  negligible. 

From  the  definition  of  the  scattered  ray  equation 


D  =  P  +  d 


u0( 


a, 6) 


(5) 


where  5  is  the  location  of  a  point  at  a  distance  d  from  the  scat¬ 
tering  point,  we  can  first  find  the  distance  dg  along  the  scat¬ 
tered  ray  at  which  the  ray  is  closest  to  the  center  of  the  earth 
as  follows. 

The  square  of  the  distance  from  the  center  of  the  earth  to 
a  location  5  on  the  scattered  ray  is  given  by 

D2(d)  *  5*5  =  P*P  +  2dP*uQ  +  d2  (6) 

The  shortest  distance  from  the  center  of  the  earth  to  a  point 
along  the  ray  5(d)  occurs  at  the  point  dg  where  the  first  deriva¬ 
tive  of  Equation  (6)  with  respect  to  d  is  equal  to  zero,  i.e., 


0 


2  P1 


V 


a, 6)  +  2d, 


(7) 


and  Dg  *  D(dg)  is  given  by 


If  Dq  <  a  and  dQ  >  0  then  the  scattered  ray  intersects 
the  surface  of  the  earth.  The  receiving  locations  where  the  rays 
(parameterized  by  a  and  6)  intersect  the  surface  of  the  earth  are 
those  that  satisfy 


a2  -  P*P  +  2  RrP‘uQ(a,8)  +  Rr2 


where  Rr  *  |P  -  R|  is  the  distance  between  the  scatterer  at  P  and 
the  specular  receive  location  at  R  .  This  distance  is  given  by 


Rr(a,B)  -  -  P*uQ  (  a,  8)  -  /a  -  DQZ(a,0) 


provided  P*un  <  0  ,  and  D2(a,8)  <  a2  . 


The  receive  locations  (a, 8)  satisfying  (10)  are 


R(a,B)  *  P  +  Rr(a,B)  u0(a,8) 


*  a(cos0  cosifr..,  cose  sin$  ,  sin@  1 
r  r  it  it  it 


where  the  geomagnetic  latitude,  0r,  and  longitude,  <frr,  are  found 
from 


tan0. 


Pz  +  Rr  uz 


(12) 


-  ,•  ,*  r  / 


tan* 


P  +  R  u 

y  -?  y 

P  +  R  u 

x  r  x 


A 

with  uQ  -  (ux,  u  ,  ufc) 


Figure  2  shows  the  locus  of  all  specular  receive  locations 
for  a  transmitter  located  at  geographic  latitude  of  61 °N  and 
longitude  150°W  (near  Anchorage,  Alaska)  and  a  few  scatterers  at 
an  altitude  of  110  km  height  (E-layer)  at  geographic  latitudes 
and  longitudes  which  are  within  the  auroral  oval  (say  at  geo¬ 
magnetic  latitudes  of  65°N,  67°N  and  69 °N).  Rays  (radiowaves) 
scattered  by  irregularities  at  geomagnetic  latitudes  greater  than 
71  °N  and  110  km  altitude  do  not  intersect  the  surface  of  the 
earth.  Similarly  rays  scattered  by  irregularities  at  220  km  or 
higher  (F-layer)  and  geomagnetic  latitudes  coinciding  with  the 
auroral  oval  do  not  intersect  the  surface  of  the  earth  either. 
If  the  transmitter  location  is  moved  to  a  location  at  a  geogra¬ 
phic  latitude  of  65°N  and  longitude  of  148°W  (near  Fairbanks, 
Alaska)  then  none  of  the  rays  scattered  by  E-layer  or  F-layer  ir¬ 
regularities  within  the  auroral  oval  intersect  the  surface  of  the 
earth.  This  implies  that  specular  scatter  from  irregularities 
within  the  auroral  oval  should  not  be  observable  at  Fairbanks, 
Alaska. 

For  comparison  purposes  we  show  in  Figure  3  the  locus  of 
all  specular  receive  locations  for  a  transmitter  located  at  geo¬ 
graphic  latitude  of  35®N  and  longitude  120*W  (geomagnetic  lati¬ 
tude  of  41 . 7 °N)  and  scatterers  at  an  altitude  of  250  km  (F-layer) 
at  the  geographic  locations  shown  and  which  correspond  to  geomag¬ 
netic  latitudes  between  46°N  and  56°N.  The  scatterers  may  be  ir¬ 
regularities  in  the  ionoization  of  the  F-layer  occurring  as  a 
result  of  a  high-altitude  nuclear  burst.  Rays  scattered  by 
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Figure  2  Locus  of  receiver  locations  (solid  lines)  where  specular  reflections 
from  indicated  E-layer  scatterer  locations  (x*s)  can  be  received  due 
to  a  transmitter  at  T 


irregularities  at  geomagnetic  latitudes  below  44°N  do  not  inter¬ 
sect  the  surface  of  the  earth. 

The  two  examples  of  Figures  2  a*d  3  indicate  that  even  at 
mid-latitudes,  both  transmitter  and  receiver  must  be  at  least  4° 
of  latitude  south  (in  geomagnetic  coordinates)  of  the  scattering 
region  in  order  for  the  specular  condition  to  be  satisfied.  This 
condition  can  be  relaxed  somewhat  when  one  takes  into  account 
that  the  scatter  from  field-aligned  irregularities  is  not  purely 
specular.  We  discuss  the  theory  of  scattering  by  anisotropic 
field-aligned  irregularities  next. 

2.2  ANISOTROPIC  SCATTERING  THEORY 

Auroral  scatter  at  VHP  is  a  phenomenon  by  which  energy  is 
scattered  by  elongated  irregularities  which  are  aligned  with  the 
earth's  magnetic  field.  The  irregularities  are  random  perturba¬ 
tions  in  the  refractive  index  of  the  medium  caused  by  random 
fluctuations  in  the  electron  density  at  heights  between  100  km 
and  500  km  above  the  earth's  surface.  Because  the  irregularities 
are  elongated,  the  random  fluctuations  in  the  refractive  index 
are  anisotropic.  Field  aligned  scatter  from  ionospheric  irregu¬ 
larities  also  occurs  at  night  at  equatorial  latitudes  and  else¬ 
where  after  a  nuclear  explosion  has  occurred  at  ionospheric 
heights.  This  latter  case  is  often  referred  to  as  the  'bomb 
mode'  of  propagation. 

The  total  received  power  can  be  expressed  as 


Pr  -  ptGtGr 


HI  • 

(4*r  v 


-°*23KA  o  d3? 

2  2 
RtRr 


(14) 


where  o  is  the  scattering  cross-section  per  unit  volume  of  the 
irregularities,  Rt  is  the  distance  from  the  transmitter  to  the 


scatterer,  Rr  is  the  distance  from  the  scatterer  to  the  receiver, 

Ka  is  the  two-way  absorption  loss  incurred  as  the  incident  and 
scattered  waves  propagate  through  the  D-region  of  the  ionosphere 
and  the  integration  is  over  the  portion  of  the  scattering  region 
intersected  by  the  transmit  and  receive  antenna  beams.  Note  that 
the  scattering  cross-section  per  unit  volume,  o,  differs  from  the 
radar  backscatter  cross-section,  by  a  factor  of  4*. 

If  we  assume  that  each  scattered  wave  undergoes  single¬ 
scattering  within  the  scattering  volume,  the  scattering  cross- 

section  of  the  field-aligned  irregularities  is  given  by 

o  «  —  k4  sin2x  An2s[k(i--A.  ),k(m,-m.  )  ,k (n,-n.  )  ]  (15) 

(2») 

where 

X  »  angle  between  incident  E  vector  and  the  direc¬ 

tion  of  scattering 

k  *  2*/ A 

An2  »  refractive  index  variance  due  to  presence  of 

field-aligned  irregularities  in  the  ionosphere 

S(*)  «  normalized  (unity  refractive  index  variance) 

three-dimensional  wavenumber  spectrum  of  the 
refractive  index  fluctuations  defined  later  in 
Equation  (19)  This  normalization  implies 

///  S(kx,ky,kz)dkxdkydkz  -  (2*)3  . 

A 

u^  ■  (i^,m^,n^)  »  direction  cosines  of  incident 

wave  vector. 

A 

Uq  *  (l2,m2>n2)  *  direction  cosines  of  scattered 

wave  vector. 


The  single-scattering  assumption  used  in  arriving  at  (15) 
is  valid  for  most  cases  of  interest  including  auroral  scatter. 
However  in  some  nuclear-disturbed  ionospheric  scatter  channels  it 
may  be  necessary  to  include  multiple  scatttering  effects.  In 
this  case  the  scattering  cross-section  is  given  by  (15)  with  the 
normalized  wavenumber  spectrum  S(*)  replaced  by  an  'effective' 
wavenumber  spectrum  which  includes  Fresnel  diffraction  effects 
and  multiple-scattering  effects  [Yang  and  Yeh,  1984],  Multiple¬ 
scattering  has  the  net  effect  of  reducing  the  forward  scattered 

A  A 

power  (u^  *  u0)  and  enlarging  the  cone  of  angles  where  signifi¬ 
cant  power  is  scattered.  The  results  presented  in  this  report 
will  be  based  on  the  single-scattering  assumption. 

The  refractive  index  variance  can  be  expressed  in  terms  of 
the  electron  density  fluctuations  by  noting  that  the  refractive 
index  in  the  ionosphere  is  given  by 


(n+an) 


where 


N 

AN 

fN 


/  f 2  f2/f2 

/,  n  i  f  an>, 

V  f 2  "  2  f2  lN  > 
1  —  I" 

f2 

mean  electron  density 


(16) 


perturbation  in  electron  density  due  to  field 
aligned  irregularities 

/bgN  plasma  frequency  of  ionosphere  (varies 
with  height) 


frequency 

e2/4w2eQm, 
mass  an 


e  is  the  electron  charge,  m  is  its 
is  the  free  space  permittivity. 
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The  peak  plasma  frequency  of  tne  lonospnere  occurs  ac  cne 
height  of  maximum  ionization  in  the  F2  layer  and  is  less  than 
15  MHz  so  that  at  VHF  f2>>fN2.*  The  refractive  index  variance  is 
then  approximately  given  by 


An2  Z  1  N  fANl2 

An  "  4  jr  • 


The  scattering  cross-section  of  the  irregularities  can  then 
be  expressed  as 


where  ST(*)  is  referred  to  as  the  two-dimensional  transverse 
wavenumber  spectrum  of  the  irregularities  and  SL(*)  is  the  longi¬ 
tudinal  wavenumber  spectrum. 

This  expression  assumes  that  the  spatial  correlation  of  the 
electron  density  fluctuations  is  separable.  One  such  example  is 


from  which  Equation  (22a)  follows.  Equation  (22b)  follows  from 
the  fact  that 


(m2-m1)2  +  (n2-n1)2  *  juj-uJ2  -  ( i2~l\ ^  (24b) 


Up  to  this  point,  the  theory  of  field-aligned  scatter  is 
very  general.  We  will  now  discuss  a  particular  way  (model)  of 
calculating  the  direction  cosines  based  on  the  use  of  geomagnetic 
coordinates  as  the  reference  coordinate  system  and  modeling  the 
earth’s  magnetic  field  as  the  dipole  field  of  a  sphere.  The 
model  could  eventually  be  refined  by  using  corrected  geomagnetic 
coordinates  and  by  including  higher  order  terms  for  the  earth's 
magnetic  field.  We  have  chosen  to  use  the  simpler  model  because 
it  is  relatively  easy  to  develop  and  it  contains  the  features 
which  are  relevant  to  the  theory  of  field-aligned  scatter.  Addi¬ 
tion  of  the  refinements  would  be  appropriate  if  and  when  a  need 
to  predict  accurately  the  path  loss  and  other  parameters  of 
field-aligned  scatter  for  specific  geographic  locations. 

2.2.1  Calculation  of  Direction  of  Magnetic  Field  Line 

The  calculation  of  the  scattering  angle  for  a  given 
geometry  requires  knowledge  of  the  orientation  of  the  elongated 
irregularities.  The  longitudinal  axis  of  the  irregularity  coin¬ 
cides  with  the  earth's  magnetic  field  lines.  Modeling  the 
earth's  magnetic  field  as  a  dipole,  as  shown  in  Figure  5,  we  can 
express  the  line  which  passes  through  a  height  r  at  latitude  0  as 


ffr.el 


cos 


constant 


(25) 


where  h=r-a  is  the  height  of  the  magnetic  field  line  (irregu¬ 
larity)  above  the  surface  of  the  earth. 

In  a  Cartesian  coordinate  system,  where  the  x-axis  is  on 
the  magnetic  equator  plane  and  the  z-axis  passes  through  the  mag¬ 
netic  north  pole,  the  constant  field  lines  can  be  expressed  as 


f(x,z)  * 


(x2+z2)3/2  * 


If  we  define  uf  as  a  unit  vector  tangential  to  the  magnetic 
field  line  at  a  given  latitude  0  and  longitude  $  *  0,  then 


*  (  3f  C  .  3f  1' i/ff3f^2  .  r3f-|2i1/2 

x  +  zJ/L[av)  +  J 


«  -  1/2  o  o  1/2 

[3xz  x  +  (2z2-  x2)z]/(x2+  z2)  (x2+  4z2) 


-1/2 


(l  +  3  sin20)  {3  cos0  sin8  x  +  (3  sin26-l)z} 


The  magnetic  field  is  parallel  to  the  earth's  surface  at  the 
equator  (6*0),  while  at  the  magnetic  pole  (0=90°),  the  magnetic 
field  line  is  normal  to  the  earth's  surface,  as  expected. 

The  unit  vector  parallel  to  the  magnetic  field  line  at 
arbitrary  longitude  ♦  is  found  by  rotating  the  above  vector  about 
the  z-axis  by  an  angle  so  that 


s 

1 
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uf  ■  (l+3  sin20)1</2{3  cos 0  sin0  cos*  x  +  3  cos©  sin©  sin*  y 

+  (3  sin20-l)  2} 

(28) 

A  A  A 

where  the  x ,  y  and  z  axes  are  as  shown  in  Figure  4.  Note  that 

A 

Uj  is  independent  of  r. 

2.2.2  Calculation  of  the  Direction  Cosines  and  Scattering  Angle 

Let  the  transmitter  be  located  at  geomagnetic  latitude  and 
longitude  and  tlie  rece*ver  at  geomagnetic  latitude  and 

longitude  (0r,*r)  both  on  the  surface  of  the  earth,  i.e.,  r  *=  a. 
Let  the  scatterer  at  geomagnetic  latitude  and  longitude  (0,*)  be 
at  a  height  h  above  the  surface  of  the  earth,  i.e.,  r  *  a  +  h. 

A 

The  unit  vector  u^  directed  from  the  transmitter  to  the 
scatterer  can  be  expressed  in  the  cartesian  coordinate  system  of 
Figure  4  as 

A  i  A 

Uj  *  -g-{(r  cos6  cos*  “  a  cos0tcos*t)x 

A 

+  (r  cos©  sin*  -  a  cos0tsin*fc)y 
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The  unit  vector  Uq  directed  from  the  scatterer  to  the  re¬ 
ceiver  is  similarly  given  by 


«  —  1  a 

uQ  ■  ^=-{(r  cos9  cos*  -  a  cos0rcos*r)x 


A 

+  (r  cos 9  sin*  -  a  cos0rsin*r)y 


A 

+  (r  sin0  -  a  sin0r)z} 


where 


2  2  2 

Rr  ■  r  +  a  -2ar  cosa^. 


cosa  *  sin0  sin0  +  cos0  cos0  cosf*  -*) 
v  r  r  r 


The  longitudinal  direction  cosine  can  be  determined  by 
noting  that 


r*i  *  (uo-ui),ui 


-1/2 


Rr~  (l+3sin^0)  {-2r  sin0+3a  cos0  sin0  cos0rcos( *r~*) 


-1/2 


+  a(3sin^0-l)sin0r }  %  Rt“A(l+3  sin^©)  {-2r  sin0 


+  3a  cosO  sin0  cos0tcos( *fc- *)  +  a(3sin2 0-l)sin0fc  }  (35) 
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which  after  some  algebra  reduces  to 


VL  *  l2~ll  *  2  sin(  *s/2)cos* 


-  -1/2  coso  coso 

(1+3  siire)  {3a  sine(-= — ^  +  — — -) 


.  .  sine.  sin0_ 
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The  scattering  angle,  *s,  can  be  found  from  its  definition 


V2  -  | u1-u2 | 2  »  4  sin2 ( 4s/2) 
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After  some  algebra  (38)  reduces  to 


2  (Rt+  Rr^  “2a2<  l-coso) 


where 
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coso  »  sinetsiner  +  cosetcosercos(*t-*r ) 
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The  tangential  component  of  the  direction  cosines  can  then 
be  calculated  from 

VT2  *  (m2-m1)2  +  (n2-n1)2  *  4  sin2 ( $g/2 )cos2 ♦ 

*  V2“VL2  (41) 

where  V  is  given  by  Equation  (39)  and  VL  is  given  by  Equation 
(37). 

2.2.3  Calculation  of  Average  Scattered  Power 

Field-aligned  scatter  is  caused  by  irregularities  which  are 
confined  to  a  range  of  latitudes  and  heights  above  the  surface  of 
the  earth.  This  is  certainly  the  case  where  auroral  scatter  is 
concerned  as  the  auroral  region  occurs  at  geomagnetic  latitudes 
between  65°  and  71*  North  in  the  night  sector  and  between  70*  and 
75"  in  the  day  sector  under  mildly  disturbed  magnetic  conditions. 
During  periods  of  high  magnetic  activity  the  aurora  moves  south¬ 
ward  and  during  periods  of  very  quiet  magnetic  activity  it  moves 
north  by  1  or  2  degrees  (Gassmann,  1973).  Field  aligned  scatter 
caused  by  a  nuclear  explosion  at  ionospheric  heights  will  on  the 
other  hand  occur  over  a  wide  geographic  area  which  may  be  equal 
to  a  significant  portion  of  the  continental  US. 

Since  the  scattering  region  is  confined  to  specific  geo¬ 
graphic  areas,  it  is  convenient  to  express  the  scattering  inte¬ 
gral  in  (14)  in  geomagnetic  coordinates  as 


gt(h,e,*)gr(h,e,*)B(h)sT(h,  e,*)sL(h,  e^jr^cose 


R*(h,e,$)  R*(h,8,$) 


where  h  -  r  -  a  is  the  height  of  the  scatterer  above  the  surface 
of  the  earth,  ST(»)  and  SL(*)  are  the  wavenumber  spectra  of  the 
irregularities  in  the  planes  transverse  and  along  the  magnetic 
field,  the  factor  B(z)  is  proportional  to  the  density  of  scat- 
terers  as  a  function  of  height  and  is  defined  as 


•(h)  -  sin2x(i2)2 


gt  and  gr  are  the  antenna  patterns,  hj  and  h2  are  the  heights  of 
the  bottom  and  top  of  the  scattering  region  (e.g..  E-layer  at 
high  latitudes  and  F-layer  at  low  and  mid  latitudes),  8j  and  &2 
are  the  geomagnetic  latitudes  of  the  southern  and  northern  bound¬ 
aries  of  the  scattering  region  (e.g.,  auroral  oval)  and  and  <|>2 
are  the  geomagnetic  longitudes  over  which  scattering  takes  place. 

The  outer  scale  of  the  irregularities  along  the  magnetic 
field  lines  is  long  compared  to  the  wavelength  so  that  the  spa¬ 
tial  correlation  of  the  electron  density  fluctuations  along  the 
magnetic  field  lines  can  be  assumed  to  be  Gaussian  and 


where  L  is  the  longitudinal  outer  scale  of  the  irregularities  and 


Kl  *  2k  sin  y  cos<|<  =  kVL(h,6,4>) 


with  VL(h,6,$)  given  by  Equation  (37).  Specular  scatter  occurs 
when  kL  >>  1  and  VL  *  0.  Note  that  we  have  normalized  SL  so  that 
/  SL(k)dk  *  2n . 

The  wavenumber  spectrum  in  the  plane  perpendicular  to  the 
irregularities  is  of  the  form 

Mkt)  *  sJ1  +  (KTLT)2rm/2 


*  cm(KTLTrm 


where  Lj  is  the  outer  scale  of  the  irregularities  in  the  plane 
transverse  to  the  magnetic  field,  m  is  wavenumber  spectrum  slope 
which  determines  the  frequency  dependence  of  the  scattered  field, 


Kt  «  2k  sin  —■  sin  ♦  »  kVT(h,9,$) 


with  VT(h,0,$)  given  by  Equation  (41). 

The  proportionality  constant  C  is  chosen  so  that 


h  L  ktst(rt)  ■  1 


i  .  e  . ,  C 


2»  L2  (m-2 ] 


,  m  >  2 
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The  density  of  scatterers  per  unit  height  is  proportional 
to  the  plasma  frequency  whose  height  variation  can  be  modeled  as 


,  ,  1/2 

fN(h)  «  f0E[l-(h-hE)Vb‘  ] 

»  hE_bE  <  h  <  hE+bE 

r  ,  v  2  2  1/2 

(49) 

fN(h)  -  f0F[l-(h-hF)VbJ] 

»  hF“bF  <  h  <  hF+bF 

where  fQE  and  fQF  are  the  critical  frequencies  of  the  E  and  F 
layers,  hE  and  hp  are  their  heights  and  bE  and  bp  are  their  half¬ 
widths.  The  critical  frequencies,  fQE  and  f0F,  exhibit  diurnal, 
seasonal,  solar  cycle  and  at  high-latitudes  magnetic  disturbance 
dependence . 

Substituting  Equations  (43)  through  (49)  in  the  scattering 
volume  integral  (42),  we  can  express  the  average  scattered  power 
as 


m-2  fnE 
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where 

CE2(m)  «  ( m- 2 ) L  LT2"m 


Cp2(m)  -  ( m- 2 ) L  LT2_m  (|^)p 
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(53) 

In  arriving  at  <  50)  —  ( 53)  we  have  assumed  that  the  critical 
frequency  of  the  scattering  layer,  fgE  and  fgF,  and  the 
normalized  electron  density  fluctuations  are  constant  within  the 
scattering  region.  We  have  also  assumed  that  all  rays  scattered 
from  irregularities  in  the  E-layer  suffer  roughly  the  same 
absorption  loss  KE  while  all  rays  scattered  from  irregularities 
in  the  F-layer  suffer  an  absorption  loss  of  Kp  dB.  These 
assumptions  could  be  relaxed  by  including  the  absorption  of  each 
individual  scattered  ray  and  the  critical  frequencies  and 
electron  density  fluctuations  in  the  integrand  of  Equation  (53). 

In  arriving  at  Equation  (53)  it  has  also  been  assumed  that 
the  transverse  outer  scale  of  the  irregularities  Lj.  is  much 
larger  than  the  wavelength  of  interest  so  that  KVTLT  >>  1.  The 
two  terms  in  (37)  are  proportional  to  the  power  scattered  by  ir¬ 
regularities  in  the  E-region  (90-150  km)  and  F-region  (150- 
500  km)  of  the  ionosphere,  respectively.  The  triple  integrals 
I(hE,bE)  and  I(hp,bp)  depend  only  on  the  geometry  and  polariza¬ 
tion  of  the  transmit  and  receive  antennas  and  can  be  evaluated 
given  the  geographic  location  of  the  terminals  and  the  scattering 
region.  Prediction  of  the  scattered  power  requires  modeling  of 
the  critical  frequencies  of  the  layers,  f qE  and  f0F,  the  absorp¬ 
tion  loss  for  E-layer  and  F-layer  scattering,  KE  and  Kp,  and  the 
structure  parameters  of  the  irregularities,  CE2  and  Cp2.  Models 
of  the  diurnal,  seasonal  sunspot  and  geographical  variations  of 
the  critical  frequencies  of  the  layers  and  the  absorption  loss 
are  readily  available.  The  structure  parameters  depend  on  the 

a 

electron  density  fluctuations,  (an/n)  ,  the  longitudinal  and 
transverse  outer  scales  of  the  irregularities,  L  and  LT  respec- 


tively,  and  the  transverse  wavenumber  spectrum  slope,  m.  The 
latter  can  be  determined  from  the  frequency  dependence  of  the 
scattered  power,  f"(m+2)#  a  model  for  the  structure  parameter 
CN2(m)  for  field-aligned  irregularities  in  the  F-region  has  been 
developed  by  Fremouw  and  Lansinger  (1981]  from  satellite  scintil¬ 
lation  data  at  latitudes  ranging  from  equatorial  to  high-latitude 
auroral  paths.  However  it  is  not  clear  whether  this  model  ap¬ 
plies  to  ground  based  terminals.  Furthermore  when  the  terminals 
are  on  the  ground  the  field-aligned  scatter  at  high  latitudes  is 
due  to  irregularities  at  E-region  heights  even  though  the  scat¬ 
tering  cross-section  (density)  of  the  irregularities  at  F-layer 
heights  is  greater.  The  reason  for  this  is  that  the  magnetic 
field  lines  are  nearly  vertical  at  high  latitudes  and  since  the 
scatter  due  to  elongated  irregularities  is  nearly  specular,  it 
was  shown  in  Section  2.1  that  the  energy  scattered  by  F-region 
irregularities  at  high  latitudes  is  scattered  up  and  not  back  to¬ 
wards  the  earth. 

From  measurement  of  the  average  scattered  power  on  a  high 
latitude  (auroral)  scatter  path  we  can  determine  the  structure 
parameter  CN2(m)  for  E-region  irregularities  provided  that  we  can 
measure  independently  the  absorption  loss  Kg  and  the  critical 
frequency  of  the  E-layer,  fgE,  in  the  scattering  region.  The 
validity  of  the  assumption  that  field  aligned  scatter  at  high 
latitudes  is  due  to  E-region  irregularities  could  possibly  be  de¬ 
termined  from  measurement  of  the  delay  power  impulse  response.  If 
both  E-region  and  F-region  scatter  occur,  then  the  delay  power 
impulse  response  should  exhibit  two  peaks  corresponding  to  the 
delay  of  the  energy  scattered  by  irregularities  near  the  maxima 
of  the  two  layers. 

The  model  of  Equations  (50)-(53)  can  also  be  used  to  pre¬ 
dict  the  strength  of  field-aligned  scatter  due  to  ionization  ir¬ 
regularities  caused  by  a  high-altitude  nuclear  explosion.  To  do 
this,  it  would  be  necessary  to  determine  the  geographical  loca- 


tion  of  the  nuclear  explosion,  the  rms  electron  density  fluctua¬ 
tions,  the  mean  electron  density  or  plasma  frequency  of  the  back¬ 
ground  ionization  after  the  explosion,  the  outer  scales  of  the 
irregularities  along  and  transverse  to  the  magnetic  field  lines, 
and  the  D-region  two-way  absorption  loss.  Computer  codes  which 
predict  these  parameters  are  described  in  the  classified 
literature . 


2.2.4  Calculation  of  RMS  Delay  Spread 

In  addition  to  the  strength  of  field-aligned  scatter,  it  is 
of  interest  to  determine  the  coherence  bandwidth  or  delay  spread 
of  the  scattered  signal.  The  rms  delay  spread  can  be  calculated 
by  noting  that  the  delay  t  of  a  scatterer  at  a  distance  Rt  and  Rr 
from  the  transmitter  and  receiver  respectively  is  given  by 


T 


V  Rr 


(54) 


where  c  is  the  propagation  velocity. 

Then  we  can  define  the  delay  moments  of  the  scattered  field 
as 


Tn  = 
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111  t-E-s — -  A 

»  c 


2  2 
RtRr 


a  exp{-0.23  K.  }d3r 

III  - 2-2^ - 

v  R3  rJ 


(55) 


where  o  is  the  scattering  cross-section  per  unit  volume  of  the 
irregularities  defined  in  previous  sections  and  the  integration 
is  over  the  scattering  volume  as  outlined  in  Section  2.2.4. 


The  rms  delay  spread  is  now  defined  as 


(56a) 


and  the  coherence  bandwidth  is  then  defined  by 


Bc  -  74-  (56b) 
3.  AURORAL  SCATTER  RESULTS 

In  this  section  we  apply  the  theoretical  model  developed  in 
Section  2  to  investigate  the  characteristics  of  VHF  anisotropic 
scattering  by  field-aligned  irregularities  in  the  auroral 
E-layer.  In  particular  we  will  determine  the  aspect  sensitivity 
of  anisotropic  scatter,  the  effects  of  polarization,  antenna 
height  and  beamwidth  and  frequency  dependence. 

In  order  to  determine  these  effects  we  need  to  establish 
the  location  of  the  auroral  E-layer  which  coincides  with  the  lo¬ 
cation  of  the  auroral  oval.  Figure  6  shows  the  location  of  the 
auroral  oval  (solid  lines)  in  (corrected)  geomagnetic  coordi¬ 
nates  under  quiet  magnetic  activity  when  the  magnetic  index 
Kp  ■  1.  The  auroral  oval  encircles  the  geomagnetic  pole  at  a 
constant  geomagnetic  latitude  during  the  day  and  moves  southward 
at  night.  The  southern  boundary  of  the  auroral  oval  also  moves 
southward  as  the  magnetic  index  increases  (1°  per  integer  in¬ 
crease  in  Kp) . 


NOTE:  The  geomagnetic  north  pole  has  been  assumed  to  be  at 

latitude  78.3°N  and  longitude  69°W. 


Figure  6  Polar  ionosphere  features  as  a  function 
of  latitude  and  local  time 


The  ionization  in  the  E-layer  in  the  auroral  region  is  due 
to  two  sources:  one  is  electromagnetic  solar  radiation  and  the 
other  is  precipitating  electrons.  The  latter  is  present  all  the 
time  and  is  the  source  of  the  so  called  auroral  E-layer. 
Electromagnetic  solar  radiation  produces  ionization  in  those 
parts  of  the  ionosphere  that  are  illuminated  by  the  sun.  Hence 
it  only  produces  ionization  in  the  day  sector.  The  critical  fre¬ 
quency  of  the  E-layer  can  then  be  expressed  in  terms  of  the  con- 
contributions  from  the  two  sources  of  ionization  as  [Gassmann, 
1973] 


f0E  =  [(f0E)4  +  (fgE)4] 


1/4 


uv 
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where  (fnE)  is  due  to  ionization  caused  by  the  solar  illumina- 
u  uv 

tion  and  is  a  function  of  geographic  loation,  time-of-day, 

season,  and  sunspot  number,  while  (f0E)  is  due  to  precipitating 

electrons  in  the  auroral  oval  and  is  a  function  of  the  magnetic 

index.  A  model  for  predicting  ( fnE)  given  the  geographic  loca- 

uv 

tion,  time,  month  and  sunspot  number  is  described  by  Malaga 
[1981],  while  Gassmann  [1973]  has  developed  a  model  for 
predicting  (fgE)  as  a  function  of  the  magnetic  activity  as  mea- 
sured  by  the  magnetic  index  Kp.  Both  of  these  models  have  been 
combined  according  to  Equation  (56)  to  predict  the  critical  fre¬ 
quency  of  the  auroral  E-layer,  fgE. 

In  the  remainder  of  this  section  we  present  the  results  of 
calculations  of  the  relative  path  loss  of  auroral  scatter  for 
different  conditions.  The  total  path  loss,  LT,  is  defined  as 

t  ._Jk _ Lr 

T  PtGtGr  C„2exp(-0 . 23  Kj 


(58) 


where  Lr  is  the  relative  path  loss  which  accounts  for  all  propa¬ 
gation  effects  except  for  absorption  loss  Kg  and  the  unknown 
structure  parameter  CE  of  the  field-aligned  irregularities. 

3.1  ASPECT  SENSITIVITY  OF  AURORAL  SCATTER 

Contour  plots  of  constant  relative  path  loss  are  shown  in 
Figure  7  as  a  function  of  geographic  coordinates  for  the  case  in 
which  the  operating  frequency  is  45  MHz  and  the  transmitter  is  at 
a  geographic  latitude  of  61°N  and  longitude  of  150°W  (Anchorage, 
Alaska).  The  calculations  assume  antennas  with  broad  azimuth 
beamwidths  and  horizontal  polarization.  The  effects  of  the  an¬ 
tenna  heights  on  the  elevation  pattern  are  included,  however. 
The  antenna  heights  assumed  are  4  meters  which  result  in  a  broad 
elevation  beamwidth  with  a  null  at  near  grazing  elevation  angles. 

The  magnetic  index,  Kp,  assumed  was  2  which  places  the 
auroral  oval  at  geomagnetic  latitudes  between  65°N  and  71°N 
during  the  night.  The  critical  frequency  of  the  E-layer  was 
f0E  *  3.87  MHz,  The  height  of  maximum  isonization  was  115  km  and 
the  bottom  and  top  of  the  layer  were  at  a  height  of  90  km  and 
140  km,  respectively.  The  field-aligned  irregularities  were  as¬ 
sumed  to  have  a  longitudinal  outer  scale  (i.e.,  outer  scale  along 
the  magnetic  field  lines)  of  10  km  which  results  in  fairly  as¬ 
pect-sensitive  scatter.  Longer  outer  scales  result  in  even  more 
aspect  sensitivity,  as  the  scatter  becomes  more  specular  the 
greater  the  ratio  of  longitudinal  outer  scale  to  wavelength  is. 
Smaller  outer  scales  result  in  more  isotropic  scatter. 

These  results  should  be  compared  with  those  of  Figure  2 
where  purely  specular  scatter  was  assumed.  The  contour  plots  of 
Figure  7  indicate  that  it  is  possible  to  observe  auroral  scatter 
in  regions  where  the  specular  condition  is  not  satisfied,  al¬ 
though  these  signals  are  much  weaker  than  in  those  regions  where 
the  specular  condition  is  satisfied  for  a  large  number  of  scat- 


terers.  The  contour  plots  also  show  that  the  areas  where  auroral 
scatter  should  be  observed  is  still  limited  geographically  due  to 
the  anisotropy  of  the  field-aligned  irregularities.  The  size  of 
the  area  where  auroral  scatter  is  observable  depends  on  the  de¬ 
gree  of  anisotropy  (longitudinal  outer  scale).  Note  that 
stronger  signals  are  received  at  locations  due  west  of  the  trans¬ 
mitter.  However  scatter  signals  can  still  be  observed  on  north- 
south  paths  provided  both  terminals  are  south  of  the  auroral  re¬ 
gion  . 

Figure  8  shows  contour  plots  of  the  delay  spread  of  the  re¬ 
ceived  scattered  field  in  microseconds.  The  delay  spread,  which 
is  a  measure  of  the  coherence  bandwidth,  increases  as  the  volume 
of  'significant'  scatterers  increases,  i.e.,  scatterers  which  re¬ 
sult  in  near  specular  scatter  for  fixed  transmit  and  recieve  lo¬ 
cations  . 

3.2  POLARIZATION  CONSIDERATIONS 

The  strength  of  the  auroral  scatter  signal  depends  on  the 
polarization  on  two  accounts:  i)  the  sin2x  factor  in  the  volume 
integral  of  Equation  (53)  where  x  is  the  angle  between  the  polar¬ 
ization  vector  of  the  wave  incident  on  a  scatterer  and  the  direc¬ 
tion  of  the  scattered  wave,  and  ii)  the  effects  of  ground  reflec¬ 
tions  on  the  elevation  pattern  for  different  antenna  polariza¬ 
tions. 

In  Appendix  A,  we  derive  the  calculation  of  the  polariza¬ 
tion  factor  sin2x  for  different  combinations  of  transmit  and  re¬ 
ceive  polarization.  We  also  derive  the  calculation  of  the  depo¬ 
larization  of  the  scattered  field  and  the  correlation  between  the 
co-polarized  (same  received  polarization  as  the  transmitted  sig¬ 
nal)  and  de-polarized  component  of  the  received  scattered  signal. 

The  effects  of  ground  reflections  on  the  elevation  pattern 
of  the  transmit  and  receive  antennas  are  accounted  by  a  factor  of 
the  form 
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gH,V(hA'A)  *  I1  +  RH>v(^)exp{j4irhA  sin(  A)/A}  j  2  (59) 

where  hA  is  the  height  of  the  antenna,  t  is  the  elevation  take¬ 
off  angle  of  the  radiated  (incident)  wave  or  the  elevation  angle- 
of-arrival  of  the  received  (scattered)  wave,  and  RHfV  is  the 
ground  reflection  coefficient  for  horizontal  (RH)  or  vertical 
(Rv)  polarization  which  depend  on  the  electrical  properties  of 
the  ground  and  the  frequency  or  wavelength  1  of  the  wave.  The 
ground  reflection  coefficient  is  given  by 


sin ( A)  -  2  v 
RH,V(A)  "  sin(A)  +  2h|v 


(60) 
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where  eg  and  og  are  the  permittivity  and  conductivity  of  the 
ground . 

The  polarization  theory  developed  in  Appendix  A  and  the 
elevation  pattern  of  Equation  (59)  have  been  used  to  calculate 
the  relative  path  loss  for  horizontal  transmit  and  receive  polar¬ 
ization,  LjjH,  vertical  transmit  and  receive  polarization,  Lyy, 
and  for  horizontal  transmit  polarization  and  vertical  receive 
polarization,  LyH.  The  latter  is  referred  to  as  the  depolarized 
component.  Table  1  shows  the  calculated  relative  path  loss  at  a 
frequency  of  45  MHz  for  a  transmitter  at  61®N  and  150°W  (Anchor¬ 
age,  Alaska)  and  a  few  different  receive  locations.  The  calcula¬ 
tions  assume  equal  antenna  gains  and  azimuth  beamwidths  so  that 
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any  differences  in  the  path  loss  are  strictly  due  to  the  polar¬ 
ization  factor  sin2x  in  the  scattering  cross-section  of  the  ir¬ 
regularities  and  ground  reflections.  A  poorly  conducting  ground 
is  assumed,  i.e.,  eg  »  10  and  =  0.001  mhos/m.  Table  1  also 
gives  the  delay  spread  for  each  case. 

Inspection  of  Table  1  shows  that  vertical  transmit  and  re¬ 
ceive  polarization  results  in  a  smaller  path  loss  by  a  couple  of 
dB.  The  reason  for  this  is  that  at  high  latitudes  the  scatterers 
which  contribute  to  most  of  the  received  energy  are  far  off  the 
great  circle  plane  because  the  magnetic  fields  are  nearly  verti¬ 
cal.  The  polarization  factor  sin2x  has  a  null  on  the  horizontal 
plane  when  the  incident  wave  is  horizontally  polarized  and  in  the 
vertical  plane  for  a  vertically  polarized  incident  wave.  On  the 
other  hand,  horizontal  polarization  would  be  more  effective  at 
equatorial  latitudes  where  most  of  the  contributors  are  near  the 
great  circle  plane.  Table  1  also  shows  that  the  depolarized 
component  is  30  to  50  dB  weaker  than  the  co-polarized  component 
which  implies  that  polarization  diversity  would  not  be  an 
effective  way  to  combat  fading  if  one  were  to  use  field-aligned 
scatter  as  a  VHF  communications  medium. 

3.3  ANTENNA  BEAMWIDTH  EFFECTS 

VHF  Meteor  Burst  Communications  systems  employ  directive 
antennas  in  order  to  achieve  a  higher  effective  radiated  power. 
Directivity  is  achieved  mainly  in  the  azimuthal  plane.  Typical 
antennas  employed  are  4-element  yagis  which  have  3  dB  beamwidths 
of  46  degrees  IH.E.  Green,  1985],  When  the  elements  of  the  yagi 
are  horizontal,  most  of  the  directivity  is  in  the  azimuth  plane. 

Since  the  field-aligned  irregularities  which  cause  auroral 
scatter  are  present  over  a  limited  geographical  area  and  at  spe¬ 
cific  heights,  the  beamwidth  and  direction  of  the  transmit  and 
receive  antennas  will  have  a  significant  effect  on  strength  of 
the  scattered  signal. 


Table  2  gives  the  relative  path  loss  and  rms  delay  spread 
at  45  MHz  for  various  combinations  of  transmit  and  receive  anten¬ 
na  pointing  angles  relative  to  the  great  circle  plane  between 
transmitter  and  receiver.  The  transmitter  is  at  a  latitude  of 
61°N  and  longitude  150°W  (Anchorage,  Alaska)  and  the  receiver  is 
at  a  latitude  of  61°N  and  longitude  of  162®W  (Bethel,  Alaska). 
The  great  circle  distance  between  these  sites  is  around  646  km. 
The  antennas  are  assumed  to  have  3  dB  azimuth  beamwidths  of  45° 
and  only  pointing  directions  north  of  the  great  circle  plane  and 
towards  the  other  terminal  as  shown  in  Figure  9  have  been  consi¬ 
dered.  The  antennas  have  also  been  assumed  to  have  omnidirec¬ 
tional  elevation  patterns  in  free  space.  However,  the  results 
include  ground  reflections  which  introduce  lobing  in  the  eleva¬ 
tion  pattern  and  a  null  for  near  grazing  angles.  The  relative 
path  loss  and  rms  delay  spread  for  the  case  in  which  antennas 
with  omnidirectional  azimuth  patterns  is  given  in  Table  2  for 
reference  purposes  also  as  well  as  the  path  loss  when  ground 
reflections  are  ignored  (NGR) . 

When  the  antennas  have  azimuth  beamwidths  of  45°,  and  they 
are  properly  oriented,  the  increase  in  path  loss  relative  to  an 
antenna  with  omnidirectional  azimuth  pattern  is  1  dB.  Ignoring 
the  effects  of  ground  reflections  results  in  17  dB  reduction  in 
path  loss. 


3.4  FREQUENCY  DEPENDENCE 


The  theory  of  field-aligned  scatter  developed  in  Section  2 
predicts  a  path  loss  with  frequency  dependence  fm+2  where  m  is 
the  power-law  dependence  of  the  transverse  wavenumber  spectrum  of 
the  irregularities,  all  other  things  being  equal  such  as  antenna 
beamwidth,  antenna  height  in  wavelengths  and  longitudinal  outer 
scale  of  the  irregularities  in  wavelengths.*  For  a  power  law  de¬ 
pendence  of  hi  s  4  and  fixed  antenna  gains,  the  path  loss  should 
increase  18  dB  every  time  the  frequency  is  doubled.  However,  for 
a  fixed  longitudinal  outer  scale  of  the  irregularities  the 
scatter  is  more  specular  at  higher  frequencies  resulting  in 
greater  frequency  dependence  than  that  predicted  by  the  power  law 
fm+2.  The  effect  of  the  anisotropy  and  resulting  specularity  on 
the  frequency  dependence  of  field-aligned  scatter  is  shown  in 
Table  3  where  we  give  the  relative  path  loss  and  the  delay  spread 
at  30,  45,  60  and  90  MHz  for  the  same  646  km  path  of  Section  3.3. 
The  antenna  height  assumed  was  0.6A  for  all  frequencies  so  that 
the  elevation  pattern  lobing  effects  is  similar  at  all  fre¬ 
quencies.  The  deviation  in  path  loss  from  the  predicted  18  dB 
per  octave  for  a  power  law  m  *  4  is  also  listed  in  the  table. 

4.  SUMMARY  AND  CONCLUSIONS 

A  theory  for  calculating  and  predicting  the  characteristics 
of  VHF  signals  scattered  by  ionization  irregularities  in  the 
ionosphere  has  been  developed.  The  theory  is  applicable  for  pre¬ 
diction  of  the  effects  of  VHF  auroral  scatter  as  well  as  the  VHF 
'bomb  mode'  of  propagation  which  is  known  to  occur  after  a  high- 
altitude  (100-500  km)  nuclear  explosion.  In  order  for  the  theory 


*  NOTE:  The  fm+2  frequency  dependence  also  assumes  that  the 
ratio  of  the  transverse  outer,  LT,  to  the  wavelength 
is  much  larger  than  1. 
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to  be  useful  as  a  prediction  tool,  data  concerning  the  electron 
density  fluctuations  and  the  outer  scales  of  the  irregularities 
along  and  transverse  to  the  earth's  magnetic  field  is  needed.  In 
order  to  predict  the  rapidity  of  the  fading,  data  regarding  the 
velocity  at  which  the  irregularities  move  across  and  transverse 
to  the  earth's  magnetic  field  lines  is  required. 

Some  of  these  data  can  be  obtained  from  measurements  of  the 
strength  of  field-aligned  scatter  signals  as  well  as  their  fading 
rate  and  multipath  spread  provided  the  critical  frequency  of  the 
scattering  layer  where  the  irregularities  are  present  and  the  ab¬ 
sorption  loss  in  the  D-region  of  the  ionosphere  are  measured 
simultaneously . 

The  theoretical  model  has  been  applied  to  investigate 
auroral  scatter  at  VHF.  In  particular  the  results  obtained 
assuming  that  field  aligned  irregularities  occur  at  E-layer 
heights  in  a  band  of  latitudes  coinciding  with  the  auroral  oval 
predict  the  following: 

i)  Auroral  scatter  is  aspect  sensitive  with  stronger 
signals  occurring  on  east-west  paths  although,  it  is 
also  possible  to  observe  auroral  scatter  on  north- 
south  paths  provided  the  terminals  are  both  south  of 
the  scattering  region. 

ii)  Delay  spreads  as  large  as  one  millisecond  are 

predicted  on  east-west  paths  provided  the  irregular¬ 
ities  occupy  a  region  extending  from  65#N  to  71°N  in 
geomagnetic  coordinates  and  the  antenna  patterns  are 
sufficiently  broad  to  illuminate  most  of  the  scat¬ 
tering  region. 


ill)  Vertically  polarized  antennas  result  in  large 
scattered  signal  strengths  at  high  latitudes  than 
horizontally  polarized  antennas  (by  2-4  dB)  because 
the  main  scatterers  are  well  off  the  great  circle 
plane . 

iv)  The  depolarized  component  of  the  scattered  signal  is 

30-40  dB  weaker  than  the  co-polarized  component  so 

that  polarization  diversity  is  of  no  use  to  combat 
fading . 

v)  The  strength  of  field  aligned  scatter  decreases  with 

frequency  faster  than  the  power  law  dependence  f“®. 

This  is  due  to  the  fact  that  the  scatter  is  more 

specular  at  higher  frequencies  so  that  fewer  scat¬ 
terers  contribute  to  the  total  received  signal. 
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APPENDIX  A 

POLARIZATION  EFFECTS  IN  FIELD-ALIGNED  SCATTER 


A. 1  INTRODUCTION 

In  this  appendix  we  develop  the  theory  necessary  to  inves¬ 
tigate  the  role  played  by  the  polarization  of  radio  waves  in  the 
scattering  from  f ield-aligned  irregularities.  In  particular  we 
would  like  to  determine  whether  horizontal  or  vertical  polariza¬ 
tion  is  more  suitable  for  field-aligned  scatter,  how  much  depol¬ 
arization  there  is,  and  if  the  depolarized  component  is  of  the 
same  order  of  magnitude  as  the  co-polarized  component  whether  the 
co-polar  and  depolarized  component  fade  independently  or  not. 

A. 2  DEPOLARIZATION  THEORY  FOR  FIELD-ALIGNED  SCATTER 

In  order  to  determine  the  role  played  by  polarization  in 
the  scattering  from  field-aligned  irregularities  we  start  from 
the  basic  equation  for  the  field  scattered  by  a  dielectric  parti¬ 
cle  (perturbation). 

Field-aligned  irregularities  can  be  modeled  as  perturba¬ 
tions  of  the  refractive  index  of  the  ionosphere.  Thus  if  n0(z) 
is  the  refractive  index  of  the  background  ionosphere  and  tn('r)  is 
the  perturbation  of  the  refractive  index  due  to  the  field-aligned 
irregularities,  then  the  scattered  field  is  given  by  [Ishimaru, 
1978] 
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where  E(r)  is  the  total  electric  field  inside  the  scattering  vol- 
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ume  V  and  uQ  is  a  unit  vector  in  the  direction  of  the  scattered 
field,  i.e. , 
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Equation  (A.l)  is  an  exact  solution  of  Maxwell's  equations 
in  terras  of  the  field  E("r')  inside  the  scattering  volume  and  the 
perturbation  of  the  refractive  index  An(7")  .  The  exact  field 
inside  the  scattering  volume  is  not  known  and  the  perturbation  of 
the  refractive  is  a  random  function  of  the  location  within  the 
scattering  volume.  An  approximate  solution  can  be  obtained  if  we 
replace  the  field  inside  the  scattering  volume  by  the  incident 
field  E^(r')  and  assume  uncorrelated  scatter.  The  latter  assump¬ 
tion  implies  that  the  field  scattered  by  irregularities  spaced  by 
more  than  the  'correlation  distance'  of  the  perturbation  of  the 
refractive  index  is  uncorrelated.  Therefore  we  can  divide  the 
scattering  volume  into  sub-volumes  of  dimensions  of  the  order  of 
the  'correlation  distance'  of  the  refractive  index  perturbation. 
The  total  average  scattered  power  would  then  be  equal  to  the  sum 
of  the  average  power  scattered  by  irregularities  occupying  each 
sub-volume.  This  can  be  expressed  mathematically  as 
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where  the  incident  field  due  to  a  source  at  i*q  with  polarization 


vector  e^  is  given  by 
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Note  that  (A. 3)  reduces  to  (A.l)  in  the  limit  as  AV  0. 

m 

We  have  chosen  to  express  the  integral  (A.l)  as  a  discrete  sum  of 
integrals  to  facilitate  the  derivation  and  interpretation  of  the 
results  that  follow. 

When  the  source  and  observation  points  are  in  the  far  field 
of  the  scatterer  at  7  we  approximate  the  wave  functions  as 
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where  rmg  are  the  coordinates  of  the  center  of  the  mcn  sub¬ 
volume  . 
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Substitution  of  (A.5)-(A.8)  in  Equation  (A. 3)  leads  to 
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where 
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The  vector  e  is  the  polarization  of  the  wave  scattered  by 
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the  mcn  irregularity  and  is  a  function  of  the  incident  wave 
polarization  and  the  scattered  wave  direction  of  propagation. 
The  average  received  power  is  proportional  to  the  received  field 
intensity  defined  as 
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which  by  virtue  of  the  uncorrelated  scatter  assumption 
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where  S(K)  is  the  three-dimensional  wavenumber  spectrum  of  the 
irregularities,  i.e., 
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In  the  limit  as  AV  ♦  0  ,  Equation  (A.  13)  reduces  to  the 

m 

solution  for  the  average  power  scattered  by  a  random  medium  pre¬ 
sented  in  Section  2  of  this  report. 

To  determine  depolarization  effects  we  go  back  to  Equation 
(A. 9)  and  express  the  polarization  of  the  wave  incident  and  scat¬ 
tered  by  each  irregularity  as  the  sum  of  two  orthogonal  polari¬ 
zations  denoted  as  horizontal  and  vertical  polarization  respec¬ 
tively.  We  also  assume  the  axis  of  symmetry  of  the  field-aligned 
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The  equation  for  the  average  scattered  power  given  in 
Section  2,  Equation  (14)  includes  the  the  effects  of 
the  antenna  patterns  and  absorption  loss  which  have 
been  omitted  in  this  appendix  to  simplify  the 
discussion . 
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irregularities  (magnetic  field  line)  to  be  the  z-axis  as  shown  in 
Figure  A.l.  Then  the  horizontal  polarization  component  is 
defined  as  the  field  component  on  the  x-y  plane  which  is  normal 
to  the  axis  of  symmetry.  Note  that  when  the  scattering  region  is 
near  the  magnetic  pole  this  plane  is  nearly  parallel  to  the  sur¬ 
face  of  the  earth  so  that  the  E^  field  component  nearly  coincides 
with  the  horizontal  polarization  as  it  would  be  defined  on  a 
coordinate  system  whose  reference  point  is  the  transmitter  or 
receiver  location.  However  this  is  no  longer  true  when  the  scat¬ 
tering  region  is  near  the  magnetic  equator.  In  this  case  a 
coordinate  transformation  must  be  done  to  convert  from  polariza¬ 
tions  in  the  scattering  volume  to  polarizations  in  the  transmit¬ 
ter  or  receiver  coordinate  system. 

In  the  scattering  volume  coordinate  system  of  Figure  A.l, 
the  vertical  and  horizontal  components  of  the  scattered  field  can 
be  expressed  in  terms  of  the  vertical  and  horizontal  components 
of  the  incident  field  as  * 
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where  the  f™j  are  complex  scattering  amplitudes  given  by 
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In  this  appendix  we  use  suffixes  9  and  $  to  denote  the 
parallel  and  perpendicular  components  of  the  incident 
field  polarization  rather  than  geomagnetic  latitude 
and  longitude. 
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The  polar  angles  and  ^®sm'*sm^  for  the  incident 

and  scattered  waves,  respectively,  are  defined  in  Figure  A.l. 

The  horizontal  and  vertical  components  of  the  total  scat¬ 
tered  field  in  the  coordinates  of  the  receive  location  can  then 
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be  expressed  in  terms  of  the  horizontal  and  vertical  components 
of  the  radiated  field  as 
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and  the  angles  8t  and  Br,  which  describe  the  coordinate  transfor¬ 
mations  for  the  transmitter  and  receiver  respectively,  are  to  be 
determined.  At  high  latitudes  where  the  magnetic  field  lines  are 
nearly  perpendicular,  8fc  -  0  and  8r  -  0. 

The  three  matrices  in  (A. 19)  describe  depolarization  due  to 
anisotropic  field-aligned  irregularities.  When  the  irregulari¬ 
ties  are  isotropic  the  total  received  field  results  mostly  from 
contributions  due  to  scatterers  which  lie  on  the  great  circle 
plane  between  transmitter  and  receiver,  i.e.,  scatterers  for 
which  ♦im“4STn  *  *  and  hence  °21  *  ®  (from  A.  18).  Also  since 
the  scatterers  have  no  axis  of  symmetry  the  coordinate  systems 
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can  be  chosen  so  that  $t=  fJr  *  0  and  we  get  the  result  that  iso¬ 
tropic  scatterers  (on  the  great  circle  plane)  do  not  depolarize 
the  incident  wave.  The  depolarization  of  the  waves  scattered  by 
isotropic  irregularities  off  the  great  circle  plane  is  small  and 
can  be  ignored  in  practical  applications.  Also,  the  scattering 
amplitude  <  a™2  when  the  scatterers  are  on  the  great  circle 
plane  so  that  horizontal  polarization  will  in  general  result  in  a 
slightly  stronger  signal  than  vertical  polarization. 

However  when  the  irregularities  are  anisotropic,  the 
strongest  contributors  to  the  total  received  field  are  not  neces¬ 
sarily  on  the  great-circle  plane  so  that  depolarization  effects 
can  be  significant.  The  fact  that  the  received  field  is  depolar¬ 
ized  does  not,  however,  imply  that  the  two  orthogonal  polariza¬ 
tions  fade  independently  as  will  be  seen  from  the  following  dis¬ 
cussions.  Also  the  scattering  amplitudes  for  horizontal  polari¬ 
zation  are  not  necessarily  greater  than  for  vertical  polarization 
since  the  main  contributors  are  not  necessarily  on  the  great 
circle  plane. 


A. 3  STATISTICS  OF  THE  DEPOLARIZED  SCATTERED  FIELD 

In  order  to  determine  which  polarization  yields  a  stronger 
signal  and  whether  the  cross  polarized  components  of  the  received 
field  fade  independently  we  need  to  evaluate  the  average  received 
power  for  each  polarization  and  the  covariance  (correlation)  for 
the  two  polarizations,  i.e. ,  we  need  to  evaluate  the  averages 
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| Egv |  ,  |Eghp  and  EgyEsh  where  the  averaging  is  over  the 
statistics  of  the  refractive  index  perturbation.  These  can  be 
expressed  in  terms  of  the  incident  power  in  matrix  form  as 
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where  use  has  been  made  of  the  uncorrelated  scatter  assumption. 
The  function  osm  is  the  scattering  cross  section  per  unit  volume 
of  the  refractive  index  fluctuations  to  be  defined  below  and  the 
3x2  scattering  matrix  Gm  is  real  with  elements  given  by  (from 
Equation  A. 19) 
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The  scattering  cross-section  per  unit  volume  is  defined  as 


(A. 24) 


where  Sfku  ]  is  the  3-dimensional  wavenumber  spectrum  of  the 
v  sm' 

irregularities  defined  in  Equation  (A. 14)  whose  argument  depends 
on  the  incident  and  scattered  angles.  In  the  coordinates  of 

A 

Figure  A.l  the  vector  ugm  is  given  by 
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In  the  case  of  anisotropic  irregularities  with  axis  of  sym¬ 
metry  along  the  z-axis,  we  can  rewrite  (A. 24)  as 
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where  ST(*)  and  SL(*)  are  the  normalized  transverse  and  longitu¬ 
dinal  wavenumber  spectra  defined  in  Section  2  of  this  report  and 
5 

An  is  the  mean  square  refractive  index  fluctuation.  The  argu¬ 
ments  of  these  spectra  are  given  by  (in  terms  of  the  notation 
given  here) 
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From  Equations  (A. 22)  and  (A. 23)  we  can  see  that  the  covar- 
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iance  of  the  two  orthogonal  polarizations,  E  E  .  ,  does  not 

sv  sh 

vanish  in  general  since  the  matrix  elements  and  G™2  are  not 

zero  for  all  m.  Therefore  the  fading  (statistics)  of  the  two 
polarizations  is  correlated.  The  correlation  between  the  two 
polarizations  is  given  by 
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Some  decorrelation  will  occur  whenever  different  scatterers 
within  the  scattering  volume  contribute  most  of  the  power 
received  on  each  polarization.  For  example  if  the  transmitted 
signal  is  horizontally  polarized  and  most  of  the  contributions  to 
the  received  horizontal  polarization  is  due  to  the  m  *  M  scat¬ 
tering  sub-volume  while  most  of  the  received  vertical  polariza¬ 
tion  is  due  to  the  m  *  N  scattering  sub-volume,  then  the  correla¬ 
tion  between  the  two  polarizations  is  given  by 
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where  PN  and  PM  are  the  total  power  scattered  by  the  Nth  and  Mth 
scatterers  and  G^2<<  Gi2'  G22<<  G22*  Using  Equations  <A.23e) 
and  (A.23f)  to  substitute  for  and  G22  in  (A.  29)  we  get 
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provided  PN  *  PM. 

A  sufficient  condition  for  uncorrelated  polarization  fading 
is  that  the  effective  scattering  volumes  for  the  two  polariza¬ 
tions  be  different,  a  condition  that  will  not  be  met  in  most 
cases  unless  the  antenna  patterns  for  the  polarizations  are  sig¬ 
nificantly  different,  i.e.,  they  illuminate  different  parts  of 
the  scattering  volume. 

A. 4  SUMMARY 

A  theory  for  analyzing  polarization  effects  on  field 
aligned  scatter  has  been  developed.  The  scattering  matrix  G 
defined  in  Equations  (A. 22)  and  (A. 23)  can  be  used  along  with  the 
theory  of  Section  2  of  this  report  to  calculate  the  effects  of 
polarization.  The  effects  of  different  antenna  patterns  for  each 
polarization  can  be  analyzed  by  multiplying  the  elements  of  the 
scattering  matrix  G  with  the  appropriate  antenna  patterns. 
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